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Abstract 

We consider extensions of the standard model based on open strings ending on D-branes, in which 
gauge bosons and their associated gauginos exist as strings attached to stacks of D-branes, and 
chiral matter exists as strings stretching between intersecting D-branes. Under the assumptions 
that the fundamental string scale is in the TeV range and the theory is weakly coupled, we study 
models of supersymmetry for which signals of annihilating neutralino dark matter are observable. 
In particular, we construct a model with a supersymmetric R-symmetry violating (but R-parity 
conserving) effective Lagrangian that allows for the s-wave annihilation of neutralinos, once gaugi- 
nos acquire mass through an unspecified mechanism. The model yields bino-like neutralinos (with 
the measured relic abundance) that annihilate to a 77 final state with a substantial branching 
fraction (~ 10%) that is orders of magnitude larger than in the minimal supersymmetric standard 
model. A very bright gamma-ray spectral line could be observed by gamma-ray telescopes. 
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Superstring theory is a promising candidate to explain the underlying symmetries of 
nature, e.g., the probable existence and breaking of supersymmetry (SUSY). In particular, 
TeV-scale superstring theory provides a brane-world description of the standard model, 
which is localized on hyperplanes extending in p + 3 spatial dimensions, the so-called D- 
branes. Gauge interactions emerge as excitations of open strings with endpoints attached 
on the D-branes, whereas gravitational interactions are described by closed strings that can 
propagate in all nine spatial dimensions of string theory (these comprise parallel dimensions 
extended along the [p + 3)-branes and transverse dimensions). The apparent weakness 
of gravity at energies below a few TeV can then be understood as a consequence of the 
gravitational force "leaking" into the transverse compact dimensions of spacetime. This is 
possible only if the intrinsic scale of string excitations is also of order a few TeV. Should 
nature be so cooperative, one would expect to see a few string states produced at the LHC, 
most distinctly manifest in the dijet [l[ and 7+jet 0] spectra resulting from their decay. 

An attractive feature of broken SUSY is that with R-parity conservation the lightest 
super symmetric particle (LSP) is a possible candidate for cold dark matter [3|. Requiring the 
relic abundance to conform with cosmological dark matter measurements serves to constrain 
the underlying theory. A consequence is that it may be possible to detect the annihilation 
products of such particles, such as gamma rays, charged leptons, and neutrinos. 

In this Letter, we propose new processes, based in brane-world string theory, for the 
efficient annihilation of neutralino LSP's (x*^'s) into monochromatic gamma rays, Z-bosons, 
charged Ws and pairs of gluons (via ~^ 77) 1^: ZZ, W~^W~, and gg). By requiring 
that the total annihilation rate generate the measured dark matter abundance we con- 
strain the parameters of the model: the string scale, the neutralino mass, the string coupling 
constant, and two unknown dimensionless parameters that depend on the details of com- 
pactification. We then calculate the gamma-ray spectrum from neutralino annihilation in 
the central region of the Milky Way, and explore the prospects for discovery with present 
and future gamma-ray telescopes. 

The basic unit of gauge invariance for D-brane constructions is a U{1) field, so that a 
stack of N identical D-branes generate a U{N) theory with the associated U{N) gauge group. 
Gauge bosons and associated gauginos (in a supersymmetric theory) arise from strings ter- 
minating on one stack. For simplicity we consider a model with 3 stacks corresponding to 
gauge groups U{3) x f/(2) x f/(l), labeled stacks a, b, and c, respectively. 

We consider the introduction of new operators, based on superstring theory, that avoid 
p-wave suppression by permitting s-wave annihilation into gauge bosons at an adequate rate. 
To create an s-wave, both gauginos must be in the same helicity state, either left- or right- 
handed. Such gaugino pair annihilation violates R-symmetry by Ar = ±2, and is therefore 
forbidden in supersymmetric Yang-Mills theory, at least at the perturbative level. However, 
it can appear in conjunction with a SUSY-breaking gaugino mass generation mechanism. 

In superstring theory, there are no conserved charges associated with continuous global 
symmetries; consequently even with unbroken SUSY, R-symmetry can be violated by higher- 
dimensional operators, although only at certain orders of perturbation theory. The R-charge 
deficit Ar is related to the Euler characteristic of the string worldsheet, x = 2 — 2(yf — /i, where 
g is the genus and h is the number of boundaries: |Ar| < —2%, with x ^ 0- Of course, only 
SUSY-preserving interactions are allowed in string perturbation theory. Note that if SUSY 
is broken at the string level, these additional restrictions are lifted. As an example, consider 
the disk worldsheet with g = 0, h = 1, hence x = 1; which incorporates the effects of all tree- 
level interactions, including the exchanges of virtual Regge excitations. Recall that gaugino 
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FIG. 1: Left: "Genus 3/2" worldsheet with one handle, with four vertices inserted at the boundary. 
Right: Two-loop open string worldsheet with two vertices inserted at the boundary c and two at 
c' while the third one is "empty". 

and gauge boson vertex operators are inserted at the disk boundary "attached" to the 
associated stack of D-branes. In this case, Ar = 0, so that two like-helicity gauginos cannot 
annihilate into gauge bosons and the amplitude for A^A^ annihilation into an arbitrary 
number of gauge bosons vanishes at the disk level, to all orders in a' = This can be 

confirmed by using SUSY Ward identities along the lines of Ref. j5|. For a gaugino pair to 
annihilate into gauge bosons one needs a worldsheet with Euler characteristic Y = — 1- It 
can be realized in two ways: a "genus 3/2" worldsheet with g = 1, /i = lil, 0], which is 
essentially a disk with a closed string handle depicted in the left-hand side of Fig. [T], and 
g = 0, h = 3, which is a two-loop open string worldsheet jl] depicted in the right-hand side 
of Fig. [U The case of a worldsheet with three boundaries is particularly interesting ^ . If one 
inserts two vertex operators creating gauge bosons or gauginos associated with one stack of 
D-branes, say c, at a single worldsheet boundary, two vertex operators associated with stack 
c' at another boundary, while keeping "empty" the third boundary as in the right-hand side 
of Fig. [H one obtains a non- vanishing contribution to the four-point scattering amplitude. 
The corresponding effective interaction is described by the supersymmetric F-term, 

Ant = ^g! N M;' (Tr WXf^W^){i:Y W^e'^W^')\^, + c.c, (1) 

where W are the usual chiral superfields with field strengths associated with appropriate 
gauge groups and the traces are taken in the fundamental representations. Here, gs is the 
string coupling, and is the number of D-branes attached to the empty boundary. (A total 
of six possibilities in the three-stack model under consideration.) The factor of 3 is the num- 
ber of choices of the empty boundary. If the stack c' ^ c, the above contribution yields the 
full amplitude. The factor F^^-^) = 3iVF(°'3) is the genus zero topological partition function 
on a worldsheet with h = 3 boundaries. It depends on the moduli of the compact space and 
takes into account various string configurations in six internal dimensions. The correspond- 
ing amplitudes are called "topological" because they are determined by the topology of the 
compact dimensions and, unlike standard amplitudes, they do not contain any kinematical 
singularities associated with Regge excitations. On the other hand, if c' = c, the four-point 
amplitude also receives non-topological contributions from all four vertices inserted at the 
same boundary ^} In this Letter, we are mainly interested in the former case, with all 
bino-like gauginos associated with the U{1) stack and gauge bosons associated with all three 



The amplitudes induced by the "genus 3/2" worldsheets, which involve closed strings propagating in the 
handle, can be analyzed in a similar way. These amplitudes are related to (and, in some sense, they are 
"square roots" of) the genus 2 topological amplitudes in type II string theory [9|. Although they are not 
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stacks. For our purposes, it is sufficient to focus on the effective interaction term of Eq. ([T]). 
In the case where the annihilation occurs through emission of a pair of U{1) gauge bosons, 
the interaction will include unknown non-topological and "genus 3/2" contributions. Note 
that such a term can be induced not only by string physics, but also by any R-symmetry 
violating extension of the standard model at a sufficiently low energy scale. 

The interaction term contained in Eq. ([T]), relevant to s-wave gaugino annihilation, is 

with lower case Latin and Greek letters labeling gauge and spinor indices, respectively, and 
the dimensionless coupling constant 

T = 3N g^, F^^'^K (3) 
In Eq. Fj^ denote self-dual gauge field strengths in the spinorial representation, 



= = F,.<^ : F^F^ = -F^^F,, + e,^,^F^'^FP\ (4) 



In order to compute the amplitude for the gaugino pair annihilation into two gauge bosons, 
we need the wavefunctions of all particles in addition to the interaction term ([2]). We 
consider the case of massive gauginos, with Majorana mass m, without addressing details of 
the SUSY breaking mechanism. In the center of mass frame, two gauginos moving along the 
2;-axis with three-momenta k and — k, respectively, are described by the Majorana spinors: 

Q with .-(±k) = (^"J'') , .-(±k) = , (5) 



where refer to spin up and down, respectively, and the rapidity 



rj = sinh ^ ( — ^ ] . (6) 



m 

On the other hand, the polarization vectors for the gauge bosons are 



,J(,.,) = ±™!>, (7) 



where refer to helicities, k is the momentum and q is an arbitrary reference vector. In 
Eq. (171) and below, we use the notation of Ref. 13]. Since, 

the interaction term written explicitly in Eq. ([2]) couples only to ( — ) gauge boson helicity 
configurations while its complex conjugate couples only to (++) configurations. In all, there 
are only four non- vanishing helicity amplitudes: 



strictly topological, one does not expect kinematical singularities to appear in the intermediate channels. 
The only difference from the two- loop open string worldsheets is that in the genus 3/2 case gauginos and 
gauge bosons must belong to the same stack of D-branes. 
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MiX'^X''^ ^ g'^^gl') = -M++ , M{X^2X'^^ ^ g'-'g'L') ^ , (10) 

where the notation corresponds to all particles incoming. 
The effective interaction F-term ([T]) yields 

= 2 TM;' _ k) (k2 = -k) e^^h,e^^hxTT{a''''aP^) . (11) 

Using Eqs. ([5]) and ([7]), we obtain 

M— = TM^^me'' (34) ^ (12) 

where we omitted the trivial S^^^^S^^^''^ group factor enforcing identical gauge charges of the 
annihilating gauginos as well as those of the created gauge bosons. Similarly, for the process 
with fermion helicities reversed, 

M+- =-rM;^me-'^{34)\ (13) 

The two remaining amplitudes, and are obtained by complex conjugating A4 

and A1+_, respectively. 

At this point we focus on one specific assignment of stacks to boundaries. With a choice 
of binos (hypercharge gauge bosons) as our LSP, and with the assumption of relatively small 
mixing with the other U{1) subgroups in stacks a and b, the bino is largely associated with 
the U{1) stack c. Under the preceding assumption of small mixing, we note that each photon 
(Z) vertex introduces a factor of approximately sin^^F (cos^w') if inserted at the boundary 
associated with the U{2) stack b, and cos^vk {sin 9w) if inserted at the boundary associated 
with stack c. To retain the purely topological structure of the amplitude we attach the 
second boundary to stacks a, b, c and leave the third boundary empty. 

In order to compute the annihilation rate, we need the sum of squared amplitudes, aver- 
aged over the helicities and gauge indices of initial gauginos and summed over the helicities 
and gauge indices of final gauge bosons: 

AA2 _ 1 

2 M6 ' ^^^^ 

where the Mandelstam variable, s = {ki + k2Yi and WW denotes final states including 
W~^W~, ZZ, 'jZ, or 77. Near threshold (s ^ 4m^o); the total annihilation rate into the 
three SU{2) gauge vector bosons is 



3c _o / h ^ 



where p = m^o/Mg. In a similar manner. 



8c _o / h 

av\ = 
'99 

and 
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The factor ( (which in principle can take any real value) parameterizes the uncertainty 
in the XX ~^ BB amplitude because of the aforementioned non-topological components 
in the matrix element, where all four vertices are attached to the same boundary, or due 
to scattering in the "genus 3/2" configuration. Dominance of the topological component 
corresponds to ~ +1. 

We now constrain a combination of the free parameters of the model by requiring that 
neutralinos have the measured dark matter abundance The density of neutralinos that 
survives after freezing out of thermal equlibrium in the early universe is given by 



o,K'-Q.l[-^\[^\ ' 'r , (18) 

^ V20yV80/ V3 X 10-26 cm3/sy ' ^ ^ 

where x-po is the neutralino mass divided by the freeze-out temperature, g^, is the number 
of external degrees of freedom available at the freeze-out temperature, and (crf)efr is the 
effective neutralino annihilation cross section evaluated at the freeze-out temperature. The 
desired effective annihilation rate, av\y^-^ + o'vlgg + ct'^Ibb ~ (^'^')cfr — 3 x IQ-^^ cm^/s,^ 
required to generate the measured relic density, I^cdm^^ = 0.113 ± 0.003, is obtained if 



2.8 ; V0.2/ V0.5/ V Ms J 

A sizable value for F^^'^^ is not implausible. As an example, consider the magnetized brane 
model whose partition function is given by Eqs. (5.28)-(5.30) of Ref. [8|. Crudely replacing 
the sum over discrete lattice momenta with integrals, one finds that F^'^'^^ is proportional a 
product of three wrapping numbers of a D9 brane around three 2-tori. This number can in 
principle be large, thus widening the available m^o-Ms parameter space. For example, if 
= 6, then for m^o = 2 TeV, M, can be probed to 4 TeV. 

These results have important implications for ongoing and future gamma ray searches for 
dark matter. In particular, neutralinos annihilating in the Milky Way halo to final states 
containing a photon (such as 77 or 7Z) lead to a very distinctive gamma-ray line which if 
sufficiently bright could provide a "smoking gun" signature of annihilating dark matter. 

If we assume little mixing with ?7(l)'s from stacks a and b, the projection onto any photon 
[Z) in the ly^W^^ final state entails a mixing angle sinOw [cos 6w)-^ For annihilation into 
the various channels we find, 

(sin^^H^ + Ccos^^vy)^ , (20) 
[cos^Ow + Csin^wY , (21) 
2 cos^ sin^ (1 - C)^ , (22) 

(23) 
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^ In addition to neutralino self-annihilation, {av)cS can potentially include the effects of coannihilation be- 
tween neutralinos and other superparticles of similar mass. We neglect contributions from these processes. 
^ This includes the range of strin g sc ales consistent with the correct weak mixing angle found in the 



f7(3) X U{2) X U{1) quiver model [11|. 
^ We note in passing that in the minimal extension of the standard model this mixing angle is fixed and 



introduces a multiplicative factor of 0.96 into the right-hand-side of Eqs. (fn)) - pT|) [12 1 



6 



(cTf I gg is given in Eq. [16] above.) The vanishing of o"f for = +1 is a reflection of the 
symmetry in Eq. ([T]) in the topological case, where the coupling is independent of the choice 
of assignment of the stacks on the boundaries. For C, = +1, these cross sections yield an 8.3% 
branching fraction to 77. The 77 fraction is much larger than is predicted by the existing 



one-loop broken SUSY calculations [13|, |lj]. For all parameter space satisfying the measured 
dark matter abundance [3], the standard annihilation rates to 77 or 7Z are typically smaller 
than about 10~^^ cm^/s. In contrast, our model predicts crv\^^ ~ 3 x 10"^''' cm^/s, which is 



more than an order of magnitude larger than the standard SUSY result. 

For neutralinos with masses above a few hundred GeV, the H.E.S.S. observations of the 
Galactic Center (GC) [l^ can be used to probe the dark matter annihilation cross section. 
The flux of gamma rays from dark matter annihilation in the GC is given by 



871 dE^ 



nlo{r) dliij) dij , (24) 

1.0. s. 



where i denotes the final state, ip is the angle observed away from the GC, dN-y/dE^ is the 
spectrum of gamma rays produced per annihilation, and n^o{r) is the number density of dark 
matter particles as a function of the distance from the GC. The integral is performed over 
the observed line-of-sight assuming a dark matter distribution which follows the Navarro- 
Frenk- White (NEW) halo profile [l6|. In Fig.[2|, the dotted curve is the gamma-ray spectrum 
corresponding to a 1 TeV neutralino with a total annihilation rate av\tot = 3 x 10~^^ cm^/s, 
that annihilates to 77 and 7Z with branching fractions of 0.1%, which is typical for a TeV 
neutralino in the the minimal supersymmetric standard model (MSSM). For significantly 
larger branching fractions to 77 or 7Z, the prospects for detection are greatly improved. The 
solid curve in Fig. [21 corresponding to C = +1, is the gamma-ray spectrum for a neutralino 
that annihilates to 77 with an 8.3% branching fraction, and does not annihilate to 7Z. 
Unlike the case of a typical neutralino, a very bright and potentially observable gamma-ray 
feature is predicted. For example, the suggestive structure at 2.5 TeV in H.E.S.S. data from 
2004 1^ can be easily accommodated within this model. If an experiment were to detect a 



strong gamma-ray line without a corresponding continuum signal from the cascades of other 
annihilation products, it could indicate the presence of a low string scale. 

In summary, within the context of D-brane TeV-scale string compactifications, we 
constructed a model that generates a supersymmetric R-symmetry violating effective 
Lagrangian which allows for the s-wave annihilation of neutralinos, once gauginos acquire 
mass through an unspecified mechanism. The model allows for a neutralino relic abundance 
consistent with the measured dark matter density. The branching fractions to monochro- 
matic gamma rays is orders of magnitude larger than in the MSSM. A very bright and 
distinctive gamma-ray line that may lie within the reach of current or next-generation 
gamma-ray telescopes is predicted. A flux near the limit presently imposed by the H.E.S.S. 
data would strongly support a near purely topological origin for the R-symmetry violating 
effective Lagrangian. 

We thank Ignatios Antoniadis and Dan Feldman for useful discussions. H.G. and D.M. 
thank the Aspen Center for Physics for hospitality. This research was supported by the 
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FIG. 2: The gamma-ray spectrum from neutralino dark matter annihilating in the Galactic Center 
(within a solid angle of IQ-^ sr), assuming the NFW halo profile. The spectrum has been convolved 
with a gaussian of width ^E^/E^ =15%, the typical energy resolution of H.E.S.S. and other ground 
based gamma-ray telescopes. The solid curve corresponds to dark matter annihilation with Q = +1, 
for which the 77 final state has a branching fraction of 8.3%. The dotted curve corresponds to 
0.1% branching fractions to 77 and 7Z, typical of neutralino annihilation in the MSSM. In both 
cases, we considered a 1 TeV mass and a total annihilation cross section of 3 x IQ-^^ era? /s. The 
continuum portion of the spectrum arises from the decay products of the W and Z bosons, and 
gluons as calculated using Pythia. Also shown for comparison are the H.E.S.S. data 15( which are 
generally interpreted to be of astrophysical origin [l7| . 
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